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Abstract 

This paper is an electronic application to my set of lectures, sub- 
ject:'Formal methods in solving differential equations and constructing models 
of physical phenomena'. Addressed, mainly: postgraduates and related read- 
ers. Content: a discussion of the simple models of linear friction, the models, 
that have the mechanism that is based on radiation reaction. The interactions 
we will deal are based on equation arrays of the kind: 



d 2 q(t) 



dt 2 

d 2 u(t,x) _ 2 d 2 u(t,x) 



Q 2 q(t) + fcom.pl (t,q,Q) 



— c 



dt 2 dx 2 
Q(t) = <l(t)\u> 



4jcS(x - x ) (F 3rc (t, q, Q)) + flit, x) 



Central mathematical points: dAlembert-Kirchhoff-like formulae. Central 
physical points: phenomena of Radiation Reaction, Braking Radiation and 
Friction. 
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Introduction. 

A Harmonic Oscillator Coupled to an One-Dimensional Field. 

We will discuss a possible description of a detail of dynamical behaviour of one- 
dimensional newtonian particle. The detail we are here interested in is named 
'friction'. We will focus on only two aspects of this phenomenon: First, if a particle 
MOVES through a medium, e.g., through water, then a special force arises, the 
force which acts on the particle so that it brakes the particle's moving, "tries to 
stop" the particle, and the particle's energy decreases. At the same time, a special 
medium motion arises: or medium waves of this or that or other kind arise, or the 
medium becomes more warm, or the medium generates a light... In such cases as 
these, one associates these two phenomena, thinking of them as reciprocal ones, 
conceiving them as the result of interplay, and says the particle generates a a kind 
of radiation which brakes the particle's motion; so, one says about the braking or 
damping radiation. 

I am now trying to express this point in a language of mathematical formulae, 
— it is just the subject of this paper. 

The Newtonian equation of motion of the particle, that moves under the action 
of an external force, F exter nai(t, q, ■ ■ •), reads: 

d 2 q(t) 

M = F exte rnal(t,q, ■ ■ ■) 

where M stands for the mass of the particle. For example, the Newtonian equation 
of motion of the linear harmonic oscillator, which moves being subjected to an 
external complementary 13 force F ext! osc{t, q, • • •), is: 

This equation is often written as 
d 2 q(t) 



dt 2 

where 



-tt 2 q(t) + fcom P i{t, <?,••■) 



n 2 := k/M , f comp i(t, q, ■ ■ ■) := F ext , osc {t, q, ■ ■ -)/M . 
We will deal with the case, where 



that is, 



dq{i^j 

Fexternali,^: 9; ' ' ') ^Tu T^T ~t~ F co7n pl(t : * * *) 



d 2 q(t) dq{t) 

M ~dp- = -^—qT + compl ( ' qi ' 



In particular, if we deal with the linear harmonic oscillator, then 



dq(i^j 

fcompl(t, 3) - ' ") = ~^~Ql ^ fcampl,0\t: 1i ' ' ") 



where 

7 := lv/M , 

1 or, 'is a source of 

2 external to the oscillator, as a physical system 
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that is, 

^ = -2 7 ^-^) + / comp; , (t, g ,-..), 

It is just — 2^ v 9q gp , the term, by means of which one simulates the physical effect 
that one is used to naming "linear" friction. Now then, I want to interpret the 
appearence of this term as an effect of radiation reaction, an effect of an interaction 
of the particle with a field. So, I have to declare models of the fields and models of 
the interactions. 

In this paper I will discuss several models of one-dimensional particle coupled 
to one-dimensional scalar field. Primarily I am interested in the model described 
by the equation array 

d 2 q(t) 



dt 2 



= -n 



(q(t) - Q(t)) + f (t) 



d 2 u(t,x) d 2 u(t,x) , , /,.. . , ,\ „ . 

^ 2 - c 2 > - 4 1C S(x - x )[Q(t) - q(t)j + fi(t,x) 



Q(t) = u(t,x ) 



I am also interested in the model described by the equation array 
d 2 q{t) 



dt 2 

32 



= -n^(t)+ 7 iQ(t)+/o(t) 



d 2 u(t,x) 2 d 2 u{t,x) /dq(t)\ 

du(t, xq) 



Q(t) = ai - 



One can say that these are models of a point interaction. In this paper I will also 
discuss several modification of these models. All they are described by an equation 
array of the form 



d 2 q(t) 
dt 2 



^ q(t) + fcompl(t,q,Q) 



d 2 u(t,x) 9 d 2 u(t,x) ./ , _,.\ , . 

^ 1 = c 2 g ( J ' - A 1C 5{x - x ) [F src (t, q, Q) j + f x (t, x) 

Q(t) = <l(t)\u> 



where l(t) stands -for a functional □, for any t fixed; in this formula we consider t 
as a free variable til . After indicating d'Alembert-Kirchhoff-like formulae for 

3 47cp = k , Q 2 = k/M ; the constants of the model mean, e.g.: c = propagating waves velocity, 
p = 'a density' of the field, k = elasticity constant, M = the mass of the particle. Of course, we 
assume c > . 

4 I use a P.A.M. Dirac's "bra-ket" syntax and suppose that q and Q are usual (one-dimensional) 
functions of t : 

q = q(t) , Q = Q(t) , 

J thus, we deal with a family {i(t)}t of functionals ; we will normally suppose that every !(t) is 
linear, for any t fixed. Moreover, we will deal with the case where < l(t)\u > is of the form 

d 

< l(t)\u >:= aou(t,xo) + ai—u(t,Xo) 
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solutions to these systems I obtain insulated effective equations of motion of q(t) 
and then I briefly compare them. 

A few words about the THREE-dimensional particle: the case is very com- 
plicated, however we can formally reduce it to the case where ONE THREE- 
dimensional particle interacts with THREE ONE-dimensional scalar fields, or, ONE 
ONE-dimensional VECTOR field, e.g., 



dt 2 



-n 2 q t {t)+ 11 Q t (t) + / ,i(t) 




Qi(t) 



duj(t, x ) 
dt 



( here i = 1, 2 



3) 



It's all that I want here to say on this difficult topic... 
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1 Models of a Point Interaction of an only 
one-dimensional Oscillator with an only one- 
dimensional Scalar Field 

In this paper we fix measure units and let x be dimensionless position parameter, 
i.e., 

physical position coordinate = [ length unit ] x x + const . 
Otherwise a confusion can ocurr, in relating to the definition 



/oo 
S(x - x )f(x)dx = f(x ) ■ 
-oo 



We assume the standard foramalism, where 

dl+(x-x ) 
5(x-x ) = g- x 

and where 1 + stands for a unit step function (Heaviside function): 



mo-U : if 



if e > o , 



1.1 D'Alembert-Kirchhoff-like formulae 

Recall that standard D'Alembert-Kirchhoff formulae read: if 

§= c2 £ +/ ' u=u ^ x )> /=/m> (*>«) w 



and given initial data, u(0, •) and du j*'^ ) 



, then 

t=o 



u = u 



(t,x) = ^ j (f(T,X + c(t-T))-f(T,X-c(t-T)))dT + Uo(t,x) 

uo(t,x) = c+(x + ct) + C—(x — ct) 



i(„(0,, + cO + U (0,,-c t )) + l£j(^l) 



x+ct 



and where / stands for any Primitive function of /, i.e., 

df(t,x) 



dx 

Note that 



f(t,x) 



f(T,x + c(t-T))-f(T,x-c(t-T)) 

does not depend on what a primitive is one has chosen!!! Moreover, we need only 
/ and not / itself. 

Now, I specify /. In this paper I will take 

/ = -4 7C (5(z - x )(F src (t)^ + h(t,x) 

For such an / I conclude that 

/= -^cl+{x - x )(F src (t)) +fi(t,x), 
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and then I infer that 



-27 



l+(x + c(t - t) - x ) - l + (x - c(t - t) - X )j (^F src (T)jdT 



+ Y c l (fl(r,X + c{t-T))- h{T,X-c{t-T)))dT + U {t, 



u(t, x) 



2c Jo 

Denote now, to be more concise, 

u i(t,x) := y c j (ji(T,x + c(t-r))- fi(T,x-c(t-r))jdT + uo(t,x) 
and then rewrite the recent relation as following: 

u{t,x) = -27 J (l + (x + c(t-T)-x )-l+(x-c(t-T)-x fj(F src (TfjdT 
+ u 01 (t,x) 

Let us now analyse this expression. We have: if ct ^ ct + (x — x ) , and if 
ct 7^ ct — (x — x ) , then 

l + (x + c(t - r) - .t ) - l+(x - c(t - t) - Xq) 

CT < Ct + (x — Xq) f 1 , Ct — (x — X ) < T 

10 
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Ct + (x — Xq) < CT I , CT < Ct — (x — Xq) 

CT < Ct — \X — Xq\ 

ct — \x — Xq\ < CT < ct + \x — Xq\ 

ct + \x — Xq\ < CT 

T < t — \X — Xq\/c 

t — \x — Xo\/c < T < t + \x — Xq\/c 
t + \X — X \/c < T 



Hence, for t > 0, 



u(t,x) = 



rt-\x-x \/c 

-2 7 J o \F„ c (T)JdT-l+(t-\x-X \/c) 



+ u i(t,x) 



Finally, this point of the analysis has an interesting consequence: 

U(t, Xq) - U i(t,X ) = -27 J (Fsrc^fjdT 

and 



u(t, x) 



u(t — \x — Xo\/c, Xq) — UQi(t — \x — Xq\/c, Xq) , if < t — \x — Xq\/c 

, if t - \x - x \/c < 

+ u i(*, x) 
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1.2 Oscillator interacting with a scalar field 

Recall that a standard relation which one is used to describing one-dimensional 
harmonic oscillator subjected to an external complementary force F ext ,osc(t) is this: 



d 2 q{t) 



^ ?(*) + fcompl(t) 



dt 2 

fcompiif) = Fext,osc{t)/M , M = the mass of the oscillated particle 
Recall also, in this paper I will discuss several systems described by 



d 2 q(t) o2 



dt 2 

d 2 u(t,x) 2 d 2 u{t,x) 



^ q(t) + fcom P l(t,q,Q) 



c 



dt 2 dx 2 
Q(t) = <l(t)\u> 



4~/cS(x - x )(F src (t, q, Q)j + /i(t, x) 



where l(t) stands for a functional, for any t fixed; in this formula we consider t as 
a free variable □ . In this case I rewrite the recent d'Alembert-Kirchhoff relation as 
following: 

rt-\x~x a \l c 

u(t,x) = -27 J ^F src (T,q,Q)jdT-l + (t-\x-x a \/c) 
+ u i(t,x) 
I have now seen: given q and uqi, then, 



in order to obtain u(t,x) I need to obtain ONLY Q(t) =< l(t)\u > 



After this observation use the last formula for u(t, x) and then obtain 
Q(t) - ^(t)||-2 7 ^ ' (F src {T,q,QfjdT-l + (t-\x-x \/c) 

+ u 01 (t,x) > 

) ( we have to consider this expression as a function of t, x ) 
We restrict ourselves to the case, where 

d 

< l(t)\u >:= a u(t,x ) + ai—uit.xo) 

i.e., 

/ d\( ft-\x -x \/c 

Q(t) = [a Q + ai — J (F src (T,q,Q))dT-l + (t-\x -x \/c) 

+ U l(t,X ) 



6 thus, we deal with a family {l(t)}t of functionals ; we will normally suppose that l(t) is linear, 
for any t fixed. Moreover, we will deal with the case where < l(t)\u > is of the form 

< l(t)\u >:= aou(t,xo) + ai-^-u(t,xo) 
at 
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Hence 



Q(t) = (a + oti^J ^- 2 7 J (F src (r,q,QfjdT + u 01 (t,x )^ 



and we have obtained: 



d 2 a(t) 

= -ft 2 q{t) + f comp i{t,q,Q) 



Q(t) = (a + a!^jy-2j J (F src (T,q,Q^jdT + u i{t,x ) 
Q(t) = (a + I u{t,x) 



X — Xq 



rt-\x-X \/c 

u(t,x) = -2 7 J (F src (T))dT-l + (t-\x-x a \/c) 
+ Uo\(t, x) , if t > . 



The specific F src and f C om P i we will discuss are: 
(A) 

^sr-cO, 9: <3) = -72?(t) + 73^<?(t) , fcompl{t, q, Q) = JlQ(t) + f (t) 

(B) 

F src (r, q, Q) = 7 o(Q(r) - q(r)) , f comp i(t, q, Q) = n 2 Q(t) + f (t) 
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2 the Models 

2.1 Effective Equation of Motion of the Particle in the Case 
of (A), i.e., in the Case where 

d 

F src (T,q,Q) = -72g(r)+73^:g(r) , fcompi(t,q,Q) = 7iQ(*)+/o(*) 



In this case the formulae 
d 2 a(t) 

= -ft 2 q{t) + f compl {t,q,Q) 
Q(t) = (a + ai^j^-2j J (F src (T,q,Q)^jdT + u i(t,x ) 

Q(t) = (ao + ai^)ju(t,a;)l 

k ) x—Xn 



u(t, x) 




(^F src (T,q,Q)jdT , if < t - \x - x \/c 

, if t - \x - x \/c < < t 



+ u i(t, x) 



become 



d 2 q(t) 

dt 2 

Q(t) = 
Q(t) = 



-n 2 q{t)+ 11 Q{t) + f (t) 

(ao+ai^) ^-2 7 y ^--/2+l3^jq(T)dT + u a i(t,x ) 
(a + J^) \ u(t,x) \ 

\ ) X—Xn 



u(t,x) = 



-27 / 
Jo 



t—\x—Xo\/c 



( - 72 + 73^:) 5 (-r)dT , \{Q<t-\x-Xo\/c 
, if t- \x -x \/c < < t 



+ u i(t, x) 
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As a consequence, 

d 2 Q(t) 
dt 2 ~ 



Q(t) +Ji(a + a i^) ^ _ 2 7 y (- 72 +73^r)<?(T)dr + M i(i,a;o)^ +/o(i) 



and then there follow quite regular transformations: 
d 2 q(t) 



dt 2 



d 2 q(t) 

dt 2 



+ 7i«4(-2 7j ( (-72+73^)?(r)^ 



+ 7i 



If 7 is a constant in £ , then 

^(0 02 „, ... . /"/ . " 



dt 2 



tt 2 q(t) - 2 7 iao7y ( - 72 + 73^)g(r)(i7 
- 271017^ -72 +73^)<?(i) 
71 (a + ai^u i(i, £o) + /o(i) 



+ 



We have just obtained an effective equation of motion of the particle subject to 
the model (A), and now, let us now try to solve this equation. We restrict ourselves 
to the case where the field is initially not excited: 

u 01 (t,x ) = 0. 

So, we now deal with the case where 

d 2 q(t) 



dt 2 



= -tt 2 q(t) - 2jia a j j ( -72 +73^)9(^)^7 



27i«i7( - 72 + 73^)<?(t) + fo(t) 



and we see, in the relation written down, the term that one writes tending to express 
the idea of linear friction: it is 

-27iai77 3 — q(t) . 
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But another detail attracts attention: this relation is not an ORDINARY differential 
relation whenever 

o ^0. 

It is because of the term 



-271007 J (-72+73^)9(1")^ 



Only if oo = we see a relation which construction is habitual: 



n 2 q(t) - 2 7l a l7 ( - 72 + l3-^)q(t) + fo(t) 



In other cases we deal with the particle's motion that one is used to qualifying as 
motion ( or, evolution ) with memory. Even if 

72 = 

the equation of motion is not ORDINARY differential equation because of q(0) in 

-271 o J l3-^q(T)dT = -27ia 73 (<?(*) - <?(0)) ( if 73 = const) 

One can say we have models of dynamics with "on only one instant concentrated 
memory". For a contrast, the case where 

72 ¥= , 73 = , 

can be referred to as a case of "wide memory" 

Let us now consider some particular cases that represent (as we think) the most 
typical properties of the general case. For simplicity, we assume all o-s and 7-s to 
be positive, and constant in t. 



"Habitual" case: 
d 2 q(t) 

dt 2 



-n 2 q(t) - 2 7l a l7 ( - 7 2 + TS^)q(t) + fo{t) . 



There is no surprise, with the possible exception of the case, where 
(7 1 a 1 77 3 ) 2 - fl 2 + 2710:1772 > , 

-7i a i773 + V (7i"i773) 2 - ri 2 + 27101772 > 

In such a case we observe self-acceleration of the oscillator, a phenomenon discussed 
and been discussing at least in electrodynamics. 

The second case is the case of " on only one instant concentrated memory" : 

= -n 2 q(t)- 27iOo773(9W -q(0)) 
d 

- 27101773 — q(t) + f (t) 
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As we have recently pointed up, this equation is not ordinary differential one. 
However, the machinery of the ordinary differential equations does here quite for. 
We illustrate it by the example, where 

/o(t)=0, 

i.e., where 

d 2 q(t) 



= -n 2 q(t) - 2 7 ia 773 (q{t) - g(0)) 
d 

dt 1 



dt 2 

- 27101773 — q(t) 



We emphasise: this equation is linear homogeneous in q(t), but write it as 
—^T + 2 7i a i773^l^ + + 27io;o773)gft) = 2 7 ia 77 3 <?(0) 



and consider this equation as inhomogeneous one. 
Thus, we infer that 



«(*) = 27ia ° 773 g( ) + e -T^^^t (Cc cog(r! t) + Cs gin(0 t)) 
+ 271^0773 



where 



fig := + 27ia 773 - (7i«i773) 2 
At first, we can calculate C c . Actually, we have 



Cc = 9(0) - . 1(0) 



i.e., 

27i«o773 
n 2 + 27100773 

Thus we have calculated C c . As for C s , it can be similarly calculated. 

An interesting detail is: if t — > +00 then g(f) has a limit, in this sense q(t) 
behaves as the usual damped oscillator. But 

q(t) -> ^ O 27iao773 — q(0) + every time that q(0) ^ and 7 ia 773 + !!! 
fi 2 + 2 7 iQ!o773 

We observe an element of a plastic behaviour! 

The last particular case we wish to discuss is the case of "wide memory": 

d 2 a(t) f* 

= -Q 2 q(t) + 2 7ia0 772 y q(r)dT 

+ 2 7:L ai772<7(t) + fo(t) 



Technically, this is the most complicated case, with the possible exception of the 
general one, for that reason we restrict ourselves to case where 



fo(t) = 0, ai =0, 
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i.e., where 



d 2 q{t) 

dt 2 



fl q(t) + 27iQ!o772 / q(r)dT 



o 



and concentrate only on the system's behaviour at large t . In order to estimate 
the asymptotic behaviour of q{t) as t — ► +00, let us handle with the characteristic 
polynomials. It is: 

A 2 + ft 2 - 271^0772 \ = (A 3 + ft 2 A - 2 7l a 772)/A 
A 

The situation is dramatic. The polynomial 

A 3 + fl 2 X - 27ia 772 
has one pure real root Ai and two complex-conjugated ones: A 2 , A 3 , A 3 = A 2 . Since 
Ai + A 2 + A 3 = , AiA 2 A 3 = 27ia 772 

we have 

Ai +2i?eA 2 = 0, Ai|A 2 | 2 = 27^0772 

and hence 

Ai > (!!!) ,i?eA 2 = i?eA 3 < 0. 

Thus, we expect an EXPONENTIAL GROWTH of the oscillator amplitude, as 
t — > +00. Does this exponential growth really exist? Actually, any expression 

Cie Al * + C 2 e A2t + C 3 e^ t , 
as a function of t, satisfies the relation 

^ - -^ + 27^0772^). 

But it in itself does not mean that such an expression satisfies the proper relation 
d 2 q(t) 



dt 2 



-£l 2 q(t) + 271CX0772 / q(r)dT . 

Jo 



A priori, all that we can now assert is that if a function, F of t , is of the form 
F(t) = de Alt + C 2 e A2t + C^ 2t , 

then 

d 2 F{t) 



8t 2 



= -tt 2 F(t) + 2710:0772 / F(r)dT + a Constant . 

Jo 



Fortunately, the problem, that has just arisen, is not difficult. Of course, every 
proper q(t) is of the form 

Cie Alt + C 2 e A2t + Cge 17 * . 

Moreover, 

q (t) = Cie Alt + C 2 e A2t + e 17 * 
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because q(t) is a particle position, hence, q{t) is real. Thus, if we take into account 
that 

t 

Cie AlT + C 2 e X2T + C 2 e X2T dr 




we can conclude: if 
then 

q(t) := C 1 e Alt + C 2 e A2t +^e^*, 

satisfies the proper relation. Thus, if we take C 2 so that 

M.(g)^0, 

e.g., C2 = A 2 , then 

x -2i?e (-^ J e Alt as t +00 . 

Thus, we really observe an EXPONENTIAL GROWTH of the oscillator am- 
plitude, as t — > +00, a factor that can throw Physicist's mind into confusion. We 
cannot here hope we have simply confounded the 'time directions'. If we had, we 
would have two roots with strictly positive real parts! In any case, we cannot 
escape from phenomenon of "self-acceleration". We defer the more detailed dis- 
cussion on this subject and notice only, that a related phenomenon is known in 
electrodynimics, see Abraham-Lorentz-Dirac equations. 

Now then, we have obtained an effective equation of motion of the particle 
subject to the model (A) and discussed properties of this model, and we turn now 
to the model (B) with the same intention. 
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2.2 Effective Equation of Motion of the Particle in the Case 
of (B), i.e., in the Case where 

F src (r, q, Q) = 7o(<2(r) - g(r)) , f comp i(t, q, Q) = Q 2 Q{t) + f (t) 

In this case the formulae 
8 2 q{t) 



8t 2 

Q(t) 



= -Q q{t) + fcom P l(t,q,Q) 



Q(t) = 



( ao + ai j^_27jf (F src (r,g,Q))rfr 



+ u i(t,x ) 



u(t, x) 



ft— \x — xo\/c 



-27 



J (Fsrc(T,q,Q)^dT , H0<t-\x-x \/c 

, if t- \x -x \/c < < t 



+ u 01 (t, x) 



become 



d 2 q{t) 

dt 2 
Q(t) 



= -n 2 q(t) + n 2 Q(t) + f (t) 



= (a + ai^j(-2j J (jo{Q(r) -q(T)]jdT + u 01 {t,xo) 



Q(t) = (a + ai^)|«(t,a:)| 



u(t, x) 



r t-\x-x \/c 

I (lo(Q(r) - q(r)))dT , if < t - \x - x Q \/c 

, if t- \x -x \/c < < t 



+ u 01 (t,x) 
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Let us focus our attention firstly on 
d 2 q(t) 



dt 2 



n 2 q (t) + n 2 Q( T ) + f (t) 



Q{t) = (ao + aij^) ( -27 f (jo(Q{T)-q{T))^dT + u 01 (t,x ) 



Write it as 



d 2 q(t) 
dt 2 



-f (t) = -n 2 (q(t)-Q(r)) 



Q(t) + (a + oti^2~{ J ^7oQ(r)^dr = (a + |^ 2 7 j (loqirfj dr + u i{t, x ) 

We have now obtained an equation array for Q and q. The next step is to form 
an insulated equation for q. LJFor this purpose, we apply, at first, the operator 

d 



to the former equation, i.e., to the equation 

cPq 
dt 2 

Then we infer 



(hl ' fo = -tf(q-Q) 



(f-«W»-i) ! 'f»(^-"(')> 



>d 2 q{t) 

v dt 2 



d_ 

dt 



'd 2 q{r) 



n 2 (q(t) + (ao + a^^j J 7o<?(r)dr 
- & (Q(t) + (a + c*i ^) 2 7 j 7oQ(r)dr) 



- /o(t) 



aa+ai §ihJ 70 ( 



d 2 q(r) 
dr 2 



fo(t))dr 

-n 2 (q(t) + (a + a 1 ^2-yj 7og(r)di 



+ fl 2 ((a + ai — ^27 j 7 g(r)dr + u al (t,x )j 



( ^_ /o(t)) + (ao+ai | )27 | 7o( ^_ /o(t))dT 

= -n 2 (q(t) +u i(t,X fj 



7 as for an insulated equation for Q, one can find it in Appendix A. 
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((/ + #)(?- / ))(t) = -fi 2 (gW + %(Mo)) 



where 



K := 



(a + "17^)27 J lo-dT, q(t) :-- 



d 2 g{t) 
dt 2 



We can stop at this last equation, or, observing that 



t d 2 q(r) dT= dq(t) dq(t) 



dr* 



dt dt 



t=o 



we can stop at that: 



(1 + 2770*0 (^^-/o(*)) 



= -2770^0-^ ft q(t) + 2770^0 ()/ 



+ Q iioi (i, so) + 277 a / fo(i~)dT 
Jo 



t=o 



Some people prefer to write such an equation as following: 



(1 + 2 77o ai)^2 + 2 77o«o^ + ft 2 )g(<) 



2770*0 



dq(t) 



dt 



+ Wu 01 {t,x a ) + (1 + 2 7 7 ai)/o(t) + 2 7 7 a / / (r)dr 

*=0 Jo 



We have just obtained an effective equation of motion of the particle subject to 
the model (B), and now, as in the previous subsection, we are trying to solve the 
resulting equation. As in the previous subsection, we restrict ourselves to the case 
where the field is initially in unexcited state: 



"oi (Mo) = 0. 



In addition, for simplicity, we take 



70 := 1, a := 1, a% := 0. 



'general abstractions are evident! 
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So, we now deal with the case where 

}2 a dq (Q 



d_ 

dt 



dt 



+ /o(0 + 2 7 / fo(r)dr 



We start out emphasising that the homogeneous equation, connected to this 
equation is exactly 



/ d 2 n d 

( 1- 9.<v — 

\dt 2 



+ 2~f- + n 2 )q(t) 
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dq(t) 



dt 



t=0 



and NOT 



d 2 



d 



(w 2+2 ^ t + Q2 )^ = 1 



In the previous subsection we have already discussed the similar factor, and we 
are using the similar machinery. The difference between the equations we have just 
now written is the rank one term 2j 



dt 



t=0 



0. This detail allows us to apply 



the usual machinery of the finite rank perturbations theory. Thus, having put 



and having taken into account the reasons of the previous subsections, we can show 
that 



<?(*) 



- *-7< 



COS(fl-yt) + 7 



sin(Q-yt) 



(?(Q) 



27 dq(t) 



n 2 dt 



t=0 



+ e 



-ft 



sin(flyt) dq(t) 



Ot 



t=t) 



27 dq(t) 



2 at 



t=0 



of course, in the case where 



/o(*)=0. 



If we now concentrate on the system's behaviour at large t , and where again, 
for simplicity, fo(t) = 0, a mathematical detail calls attention. We observe: 



q(t) -» 2 7 



Sg(t) 



c9/ 



as t 



t=o 



Again, as in the previous subsection, we observe an element of plastic behaviour. 
A surprising detail in the new situation is: the limit function 



?«>(*) = 27 



dq(t) 



dt 



t=0 



is NO solution to 



dt 2 



Ot 



27 



dq(t) 



dt 



t=0 



9 and in subsection 1.1, as well! 

10 This term, as a function of t , is a FIXED function of t, in this context a constant non-zero 
function, e.g. 1, multiplied by a CONSTANT depended on q, i.e., by a fixed functional of q. Using 



the Dirac's syntax, one can write this term as \a >< b\ with \a >= 1 and < b\q >= 2-y 



9q(t) 
at 
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every time that 

2 7 ^ 



t=o 



at 

because the q = qoc is a constant, hence its derivative is zero: 



dq(t) dqoc d t dq(t) 



dt 



t=o dt dt \ dt 



t=o 



0. 



Similar phenomena, one can detect them in the electrodynamics of moving charges. 

Now then, we have obtained effective equations of motion of the particles subject 
to the model (A) and, resp., to the model (B), and discussed properties of these 
models. 

We have seen that the linear friction can actually be described as a result of 
radiation reaction. In addition we have seen a very simple model of plastic behaviour 
of dynamical system (section 2.1) and a little more complicated model displaying 
the same effects (section 2.2). But the specific properties of these two models are 
different. In the first model memory is a function of initial position, whereas in the 
second model we rather deal with a function of initial velocity. The memory effects 
in the two models have an interesting specificity: the moving particle 'keeps in its 
memory' only initial data and 'forgets' the rest ones, with the possible exception of 
the 'past immediate': really, one needs this 'past' to calculate the derivatives of the 
position, i.e., velocity and acceleration! . . .We have also seen a model of dynamical 
system with 'wide memory' and self-acceleration . . . 

Of course, there are other interesting properties of the models, which are pre- 
sented, and many interesting abstractions and generalizations are possible. Never- 
theless, it does not form the subject of this paper. 
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3 Appendix A. An insulated equation to Q(t) 

After q(t) is found, we can determine Q{t) , at least formally, by solving 



Q{t) = [ao + a 1 ^jl-2 1 ^(jo(Q(T)-q(T)) S jdT + uoi(t,x 



or 

0(0) = (ao + ai^)uoi(t, a:o) + «i ( - 2 7 ( 7o (Q(0) - g(0))) ) 

t=o V / 

Thus we have already reduced our model, a model of an oscillator coupled the a 
scalar field, to a pair of linear 'ordinary' differential equations. Nevertheless we want 
to continue to analyse the matter and we now go searching for another relationships, 
which would simplify calculations of q, Q and u. 

At first, we will obtain another insulated equation for Q, differently and in a 
different form. 

We have 

d 2 o(t) 

-^-fo(t) = -fi 2 (g(t)-Q(r)) 
Q(t)+ (a + ai-|)27^ ( 7 0<3(r))dT 

= (ao + ai^j^ J (j q(T^jdT + u 01 (t,x )j , 



i.e., 



d 2 g{t) 
Ot 2 



-fo(t) = -&(q(t)-Q(T)) 



Q(t) - u 01 (t,x ) + (a + aij^2 7 ^ 7o (q(t) - u i(r, Zo)) c?t 

= (ao + a!^)^ j Jo(q(t) -u i(T,x fjd,T 

Denote now, to be more concise, 

Qd(t) == <3(*) - u Q1 (t,x Q ) , D i := (a + cti^J 



Then we infer 



(J^ + fl 2 ) (g d (t) + An2 7 ^ 7o g d (T)dr) 

9 2 



= (op + n2 ) D oi hi 7o(g(r)-«oi(T,a:o)) 



dr 
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and 

( + f! 2 ) (Q d {t) + An 2 7 / loQd(r)dT 



Note 



d 2 r d r d 2 d 

— 2 ^ l0 q(r)dr = - fm (t) = ^ ^ q(r)dr + - lm (t) 



t=0 



Then 

-)2 



= D oii^l (^ + r» 2 )7og(r)dT + 27^ 7og (t)|^ o - (^ + f! 2 ) (27^ 70^01(^x0)^) 



Then 

-)2 



(^2 +tf)(Q d {t)+D m 2 1 J Qd(T)dT^ 

= An (27^ 7o(« 2 0(T)+/o(r))dT + 27^7o g (t)| 4=o - (^ + r! 2 )(2 7 ^ 7o«oi(r, ar )dr 



Then 

32 



(^ + n 2 )(Q d {t) + D m 2 1 J l0 Q d {r)dr) 

'01 (27^ 7o(« 2 0(r)-0 2 Uoi(t, x )JcZt + 27— Joq(t) 
+ 2 7 ^ 7o/o(r)drj 



(27 / 7 u i(t, £0)* 



Then 

-)2 



( + fl 2 )Q d (i) + + " 2 )2 7 ^ JvQd(r)dT 

= ^01(27^ 7o(^ 2 Q(r) -O 2 M i(T,x ))dT + 27^7 0g (t) 



t=o <9i 2 



( 2 7 y 7o"oi(t, x )dr^ 



+ j 70/0 (r)dr 



Use now that 



- uoi(*,a;o) =: Qd(t) 
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Then 

-)2 



(rt g Q2 . rt 

27 J ~/ a n 2 Q d (T)dT + 27— j Q q(t) ^ - [2j J j a u i(T,x a )di 

+ 2 7 ^ 7o/o(r)drj 



Then 



= An I 27— 7oQ-(i) 



t-o ~~ ^( 27 / ^° Uoi ( T ' a;o ) dT ) +27 / ^°f°( T ) dT 



Then, finally, 



Vat 2 



+ £>oi 27^70 + ^ 2 )QdW 



= An 27 



dt 



t-o ~ 7 9i J l0 ^ T > 



( d 

where Qd(t) := Q(t) - u i{t,x ) , D i := [a + ai — 



On the surface, this equation appears to be a second order ordinary differential 
equation. It is not exactly the case. We may not arbitrary take the initial data for 
Q(t) because 

Q(t) = (ao + ai A)^-2 7 ^ ( 7o (Q(r) - g(r))) dr + u i(*,*o)^ 



and 

at 



= (a + ai^) ^-2 7 ( 70 (Q(t) - g(t))) +«oi(t,a;o)j , 

g(0) = (ao + ai^)w i(t, a*) + «i ( - 2 7 ( 7o (Q(0) - g(0))) 

t=o V 



The proper ones are these: 
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dQ(t) 



dt 



(a + ( - 27(7o(QW - q(t))) + uoi(*> flfo) 



t=o 



t=0 



0(0)- (a + ai^)woi(t,aro) +aj - 2 7 ( 7o (Q(0) - g(0))) 

t=o V 



Thus, we have already obtained two simple 'ordinary' differential equations for 

q(t) , Q(t) 

If we have found these quantities, we can try to find u(t, x) by 



rt-\x-x \/c 

u ( t ,x) = J ~ 2 7 / (^F src (T,g,g)JdT , if < t- \x -x \/c 

, if t- \x -x \/c < < t 



+ u i(t,x) 
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4 Appendix B. Remark. Complete Reflection 

Suppose, we DO have a situation where Q(t) = . In that case 

(a + ^ - 27 y (F src (T,g,Q)^«iT + ?ioi(t,a;o)^ =0 

and then 

-27 J (F src (T,q,QfjdT + u 01 (t,x ) = const e^ aoit 

for suitable constants const 0, A . Thus, we have 

u(t,x) = const e- aol ^ x - xo ^^ 

-u i{t- \x-x \/c, x ) + u 01 (t,x ) 
if < t — \x — a; |/c 

and, of course, 

u(t,x) — u i(t,x) , if t — \x — x \/c < < t 

Suppose in addition, that u i{t 7 x) is of the form u + (t + x/c) , i.e., u i(t,x) is a 
" wave moving from right to left" . Then 

u(t,x) = const e- am ^ x - xo ^ 

- u+(t - \x - x \/c + x /c) +u + (t + x/c) 
if < t — \x — x \/c 

Hence 

{constoe- "' 1 ^-^-^/^ 
const e- aol( - t -\ x - xo ^^ - u+(t - x/c + 2x Q /c) + u+(t + x/c) 
u+(t + x/c) 

if 

0< t — \x — xo\/c, x < x 
< t — \x — xo\/c, x > Xq 
t - \x - x \/c < < t 

respectively. It is because 

— u+(t — \x — xq\/c + xq/c) + u+{t + x/c) = if x < xq 
— u + (t — \x — x \/c + x /c) = — u + (t — x/c + 2x /c) if x > x a . 

In this situation we can say, that the incident wave, Uo\(t,x) = u + (t + x/c), is 
completely reflected, rejected, by the oscillator. 
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Some interesting details: 

For the F srci f CO mpi we wish to discuss, the conditions 



Q = , fa = 



mean that 



d 2 



?(*) 



dt 2 



That is, 



q(t) = A s sin(fit) + A c cos(fH) 



where A s , A c are suitable constants. Therefore, the complete reflection occurs 
where 



for suitable constants A s ,A c ,Aq . Thus, we have seen: The complete reflection 
occurs only where itoi(i, xo) is trigonometric up to constoe~ aoit ; moreover, — only 
where the spectrum of the trigonometric part of Uoi(t,Xo) contains only one 
non-zero (and real) frequence; moreover, — only where this frequence coincides 

WITH THE EIGENFREQUENCE OF THE OSCILLATOR. 



11 recall, F arc (r, q,Q) = const\q(r) + const2-^pq(r) + const^Q(r) 



uqi {t, xq) — A s sin(f2i) + A c cos(fit) + Aq + constoe 



— aoit 



26 



S.A. Chorosavin 



References 

[AK] S. Albeverio and P. Kurasov, Singular Perturbations of Differential 
Operators. Solvable Scrodinger Type Operators, London Mathematical So- 
ciety: Lecture Note Series. 271, 1999. Cambridge university press 

[BF] F.A. Berezin, L.D. Faddeev , Remark on the Schrodinger equation 
with singular potential, Dokl. Akad. Nauk. SSSR, 137 (1961) 1011-1014 
(in Russian). 

[Do] W. Donoghue , On the perturbation of spectra, Comm. Pure App. Math. 
18 (1965) 559-579 

[Fog] S.R. Foguel, Finite Dimensional Perturbations In Banach Spaces, 
American Journal of Mathematics, Volume 82, Issue 2 ( Apr., 1960 ), 
260-270 

[Fr] K.O. Friedrichs, Perturbation of Spectra in Hilbert Space, American 
Mathematical Society, Providence, (1965) 

[Jack] J.D. Jackson: Classical electrodynamics. 

John Wiley & Sons, Inc. New York-London, 1962. 

see 'Radiation Reaction, Abraham Lorentz Equation, Braking Radiation' 
and all that. 

[RT79] Razumov,A.V.; Taranov,A. Ju. 

Dipolc interaction of an oscillator with a scalar field. (Russian) 
Teoret. Mat. Fiz. 38 (1979), no. 3, 355-363. MR 80b:81034 

Electronic Print: 

[DerFr] mp_arc 02-275 

Derezinski J., Fruboes R. Renormalization of the Friedrichs Hamil- 
tonian (16K, LATeX 2e) 



\protect\vrule widthOpt \protect\href {http : //www . ma . utexas . edu/mp_arc-bin/mpa?yn=02-27 
\protect\vrule widthOpt \protect\href {http : //mpej .unige . ch/mp_arc-bin/mpa?yn=02-275}{h 
\protect\vrule widthOpt \protect\href {http : //www . maia . ub . es/mp_arc-bin/mpa?yn=02-275}{ 

[Dcr] mp_arc 02-300 

Jan Derezinski Van Hove Hamiltonians — exactly solvable models of 
the infrared and ultraviolet problem. (62K, LATeX 2e) 



\protect\vrule widthOpt \protect\href {http : //www . ma . utexas . edu/mp_arc-bin/mpa?yn=02-3C 
\protect\vrule widthOpt \protect\href {http : //mpej .unige . ch/mp_arc-bin/mpa?yn=02-300}{h 
\protect\vrule widthOpt \protect\href {http : //www . maia . ub . es/mp_arc-bin/mpa?yn=02-300}{ 

[BrDeB] mp_arc 01-275 

L. Bruneau, S. DeBievre A Hamiltonian model for linear friction in 
a homogeneous medium (548K, postscript) Jul 17, 01 



An Interaction of ... . III. Linear Friction as Radiation Reaction 27 



\protect\vrule widthOpt\protect\href {http : //www . ma . utexas . edu/mp_arc-bin/mpa?yn=01-275}-l 
\protect\vrule widthOpt\protect\href {http : //mpej .unige . ch/mp_arc-bin/mpa?yn=01-275}{httT 
\protect\vrule widthOpt\protect\href {http : //www.maia.ub . es/mp_arc-bin/mpa?yn=0l-275}{ht1 

LANL E-Print 



[AMN] Paper: lphysics/0001009 



From: adolfo@lafexSul.lafex.cbpf.br (Adolfo Malbouisson) 
Date: Wed, 5 Jan 2000 19:20:08 GMT (15kb) 

Title: An Exact Approach to the Oscillator Radiation Process in an Ar- 
bitrarily Large Cavity 

Authors: N.P. Andion, A.P.C. Malbouisson and A. Mattos Neto 
Comments: 27 pages 

Subj-class: Atomic Physics; Mathematical Physics 



[BKZ] Paper: math-ph/0210051 



From: Volkcr Bach <vbach@mathematik.uni-mainz.de> 

Date: Wed, 30 Oct 2002 00:06:05 GMT (31kb) 

Title: Mathematical analysis of the photoelectric effect 
Authors: Volker Bach, Frederic Klopp and Hcribcrt Zcnk 



Comments: See also littp://www. intlpress.com/ATMP 



Subj-class: Mathematical Physics 

Journal-ref: Adv. Theor. Math. Phys. 5 (2001) 969-999 



[CI] Paper: math-ph/0206009 



From: Dariusz Chruscinski <darch@phys .uni .torun.pl> 

Date: Fri, 7 Jun 2002 06:45:01 GMT (12kb) 

Title: Resonant states and classical damping 

Authors: Dariusz Chruscinski 

Comments: 15 pages 

Subj-class: Mathematical Physics 

A HREF="http://arXiv.org/abs/math-ph/0206009" 



[C2] Paper: |math-ph/0301024 



From: Dariusz Chruscinski <darch@phys .uni .torun.pl> 

Date: Fri, 17 Jan 2003 08:11:58 GMT (12kb) 

Title: Quantum mechanics of damped systems 

Authors: D. Chruscinski 

Comments: 20 pages 

Subj-class: Mathematical Physics 

A HREF="http://arXiv.org/abs/math-ph/0301024" 



[MPB] Paper: |hep-th/9207033| 



From: physthQulb .ac.be 

Date: Fri, 10 Jul 92 16:52:54 +0200 (20kb) 

Title: On the Problem of the Uniformly Accelerated Oscillator 
Authors: S. Massar, R. Parentani, R. Brout 

Comments: 14 pages (+postscript figures attached), ULB-TH-03/92 



